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Abstract: We investigate the universe evolution from inflation to late-time acceleration in 
a unified way, using a two-component fluid constituted from extended Chaplygin gas along¬ 
side a phantom scalar field. We extract solutions for the various cosmological eras, focusing 
on the behavior of the scale factor, the various density parameters and the equation-of- 
state parameter. Furthermore, we extract and discuss bouncing solutions. Finally, we 
examine the perturbations of the model, ensuring about their stability and extracting the 
predictions for the tensor-to-scalar ratio. 
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1 Introduction 

Cosmological observations verify the accelerating expansion of the universe, indicating that 
the transition from deceleration to acceleration was realized in the recent cosmological past 
[1, 2], The reasonable explanation would be the presence of a simple cosmological constant, 
however the possible dynamical features have led to two main alternative directions. The 
first is to modify the gravitational sector acquiring a modified cosmological dynamics, as 
in f(R) gravity [3-6], in Gauss-Bonnet gravity [7, 8], in Lovelock gravity [9, 10], in Weyl 
gravity [11, 12], in Horava-Lifshitz gravity [13-17], in nonlinear massive gravity [18-21], in 
/(T) gravity [22-31], etc (for a review see [32]). The second way is to modify the universe 
content, introducing the dark energy sector, with its simpler candidates being a canonical 
(quintessence) scalar field [33-37], a phantom field [38-44], the combination of both fields 
in a unified (quintom) scenario [45-47], or more complicated models such as in K-essence 
[48], Hordenski [49] and generalized Galileon theories [50-52] (for a review on dark energy 
see [53]). Finally, since the dynamical nature of dark energy introduces the “coincidence 
problem”, namely why are the current dark energy and matter densities of the same order 
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although they evolve differently, many extensions of the above scenarios were developed, 
with an explicit interaction between the dark energy and dark matter sectors [54-63]. 

Since in the above scenarios the dark matter and dark energy belong to two different 
sectors, one could try to construct minimalistic cosmological scenarios where both of them 
are attributed to the same source. This can be performed by assuming a cosmic fluid with 
a Chaplygin gas equation of state [64-66] or its generalization [67-75]. In particular, in this 
scenario the cosmic fluid behaves as pressureless at early times, thus mimicking the matter 
epoch, and then its equation of state progressively decreases, resulting to a cosmological 
constant at late times, and thus mimicking the dark-energy domination. However, since the 
simple models are not favored by the data [76], many extensions appeared in the literature, 
such as the modified Chaplygin gas (MCG) [77-79], the modified cosmic Chaplygin gas 
(MCCG) [80, 81], and the extended Chaplygin gas (ECG) [82-86]. 

Apart from the late-time acceleration and its previous matter epoch, a successful cos¬ 
mological scenario should be able to explain the initial phase of the universe, that is the 
inflationary exponential expansion [87, 88], or alternatively the bounce realization [89]. 
Unfortunately, the majority of the mentioned models, although successful in obtaining the 
correct late-time behavior, are unable to provide such a unified description of the whole 
thermal history of the universe. Therefore, it is both interesting and worthy to try to 
construct cosmological scenarios capable of providing such a unified description. 

In the present work, we construct a model where a cosmic fluid with an extended 
Chaplygin gas equation-of-state interacts with a phantom field. This scenario has many 
advantages, and the resulting cosmology proves to be very interesting. In particular, we 
can obtain both inflationary and bouncing solutions, thus being able to describe the initial 
universe phase in both ways. Additionally, we can obtain the subsequent matter-dominated 
phase and then the late-time acceleration, which completes the universe history. We men¬ 
tion that the phantom field may lead the dark energy sector to lie either in the quintessence, 
or in the phantom regime, or even exhibit the phantom-divide crossing during the evolu¬ 
tion, which is a great advantage. Finally, the Chaplygin gas nature, provides as usual a 
unified picture of the dark energy and dark matter sectors, offering a minimalistic picture. 
This paper organized as follows. In next section we introduce our model and write general 
equations. In section 3 we obtain cosmological solutions for various situations. We find 
the early and late time solutions, also discuss about inflationary and bouncing solutions. 
In section 4 we study scalar perturbations and investigate stability of the model. Further¬ 
more we use recent observational data to fix some solutions. Finally in section 5 we give 
conclusions. 

2 The model 

In order to construct a model of a unification of the dark sector, which is moreover capable 
of describing the phantom regime and alleviating the coincidence problem we consider a 
cosmic fluid with an extended Chaplygin gas equation of state, interacting with a phantom 
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field. In particular, the total action in the scenario at hand writes as 


S = j d A Xy/—g 


+ r 2 !f V - V{<j)) 


+ S, 


ecg j 


( 2 . 1 ) 


where n 2 = 87 tG is the gravitational constant, R is the Ricci scalar, cj> is the phantom field 
with V(cj>) its potential, and S ecg is the action of the extended Chaplygin gas. Throughout 
this work, we consider a spatially-flat Priedmann-Robertson-Walker (FRW) metric of the 
form 

ds 2 = —dt 2 + a 2 (t) 8ijdx l dx\ ( 2 . 2 ) 

where x l are the co-moving spatial coordinates, t is the cosmic time, a(t) is the scale 
factor and H = a/a is the Hubble function (dots denote differentiation with respect to t). 
Concerning the phantom held, and assuming homogeneity, its energy density and pressure 
are respectively given as usual by 

P<t> = ~^4> 2 + V(4>) 

= -V (2.3) 


Concerning the extended Chaplygin gas, in order to obtain a full generality we consider an 
equation of state of the form [82, 86 ] 

Pecg = AnPecg ~ ) (2-4) 

n Pecg 


where A n , B and n are constants. Note that in the case n = 1 the above expressions 
recovers the standard modified Chaplygin gas. 

Variation of the action (2.1) with respect to the metric gives rise to the two Friedmann 
equations, namely 

3 H~ = p ecg + p^, (2-5) 

and 

2 H — Pecg T P<jj T Pecg 4“ Pfii (2-6) 


where we set 87 rG = 1 for simplicity. In the scenario at hand, p eC g incorporates both the 
dark matter and a part of the dark energy, while p$ is added to that part of the dark 
energy and completes the dark energy sector. Finally, it proves convenient to define the 
deceleration parameter 


Q = 



(2.7) 


which quantifies the accelerated nature of the expansion (q < 0 corresponds to accelera¬ 
tion). 

Concerning the evolution equations of p $ and p ecg , the total energy momentum tensor 
conservation allows for an interaction of the form 


Pecg T 3i7 (Pecg T Pecg ) — Qi 


( 2 . 8 ) 
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and 


(2.9) 


Ptj> + 3H ( Pfj) + Prf)) — Q, 

where Q is a general interaction term (Q > 0 corresponds to energy transfer from the 
extended Chaplygin gas to the phantom field and vice versa). There are many proposed 
forms of Q in the literatures, however a general consensus is that it is reasonable to be 
proportional to some energy density, since increased energy densities are expected to lead 
to increased interactions. In the following we consider 

Q = bH(p ecg + p<p), (2.10) 

with b a constant, but other forms of interaction can be straightforwardly incorporated. 


3 Background evolution 


First of all we assume that extended Chaplygin gas density is dominant, that is we assume 
Pecg P<f)i thus the phantom field is just a small correction. 


3.1 Late-time behavior 

3.1.1 The general case 

2 

In the late time cosmology with the scale factor of the form a{t ) oc f 3 ( 1+ “) one can obtain 
the following Hubble expansion parameter, 


H = 


2 

3(1 +u)t' 


(3.1) 


with the constant uj. It yields to the following deceleration parameter, 


1 3 

q = —I— co. 
1 2 2 


(3.2) 


It is obvious that —1 < uj < —1/3 gives —1 < q < 0. 
(2.10) simplified as follow, 

Q = _ — _ 

V 9(l + w)f 3 ' 


Interaction term given by the Eq. 


(3.3) 


Analyzing of this case verifies stability of our model at the late time and confirm it. 
There is also another interesting cosmological parameter so called jerk parameter, 


a 2 d 3 a 
^ a 3 dt 3 


1 + -u(l + u). 


(3.4) 


Some observational data such as q = — 1, and j = 1 [90] suggest that u = — 1. However this 
leat to a singularity at the Hubble expansion parameter and interaction term Q. Therefore, 
we can choose uj = —0.5 to find q = —0.25, and j = —0.125 in agreement with the result 
of the Ref. [91]. Therefore, we can use the late time scale factor as follow a oc f 4 / 3 . 

By using above observational data we can obtain Hubble expansion parameter in terms 
of redshift which illustrated in the Fig. 1 for n = 1, 2,3,4. It is clear that the case of n = 1 
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coincides with the observational data (for example see the Ref. [92] ). We can see that cases 
of n = 2 and n = 3 increases value of the Hubble parameter at high redshift but decreases 
value of the Hubble parameter at present (z ~ 0). However, the case of n = 4 may be agree 
with the observational data at the high redshift. It tells that our model under assumption 
Pecg P<j> may be valid at the high redshift, means the early universe. 


(») (b) 



Figure 1. Hubble parameter in terms of redshift for A = 1/3, a = 1/2. (a) n = 1, (b) n = 2, (c) 
n = 3, (d) n = 4. 

As we know the current value of the Hubble expansion parameter is a constant in the 
range of Hq = 65 — 75 km/s/Mpc. In that case, we have the following differential equations 
corresponding to p ecg and p$ respectively, 

9 

Pecg 3 Hq pecgiX ^ ecg ) “1“ ^^0 ~ (3-5) 

and, 

Prj> — 3i7o/?b(l 3" w 0) + -Hq = 0, (3-6) 

where we should choose —1.2 < oj$ < —1 in agreement with current observational data. 
Therefore, p$ can extracted easily. Hence we could decouple conserved equation to extract 
energy densities. 
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3.1.2 The case n = 1 


This is corresponding to MCG. We can obtain two possible solutions depend on the energy 
density of ECG. If p ecg > 1, then we have the following energy density, 


p — —3Hq(5 + 


4(1 + 4) 


) + Ce 


i n 


+ 


/ 6H 0 

\J 165(1 + 4) + 9 H% + 8(1 + A)e~^ 


4(1 + 4) 


(3.7) 


where C is an integration constant. In order to obtain the above solution we used the fact 
that tanh -1 (/(p)) is not well defined for p mcg > 1. 

On the other hand, for the case of p mcg <C 1 we have the following energy density, 


P = 


- 3 "”( 5 + 4 (TW + Ce ‘"' W 

y / 165(l + 4) + 917g 


40 -4 41 


- tanh 


-i 


(1 + A)\j 165(1 + 4) + 9 Hnt 


(3.8) 


where C is an integration constant. We should note that the above solution obtained for 
a = 1. In both cases we can obtain q = —1 which is expected. 


3.2 Early-time behavior 
3.2.1 The case n = 1 

In the simplest case we consider the early universe cosmology where energy density is very 
high and assume that a = 0.5 and n = 1. Under these assumptions the ECG behaves 
as barotropic fluid with linear EoS. In that case we can obtain analytical solutions of the 
equations (2.8) and (2.9). The energy density of the ECG can obtained as follow, 

-2 

, (3.9) 


Pecg — 


^(1 + A + t)< + C, 


where C\ is an integration constant and should be of order time scale at the early universe 
to give high density. Therefore, we choose C\ <C 1 and use relation (3.9) in the equation 
(2.9) to find, 


Pcfi 


\/6 b 

((64 + b + 3) (6 + 3) + 9A 2 )t + C 2 


(3.10) 


where C 2 in an integration constant and should be of order unity because of our initial 
assumption p ecg ^$> p<f>. In that case we can obtain Hubble expansion parameter as follow, 


H 


1 I 4 ~ 36 b 2 

3y (1 + 4 + !) 2 * + (((64 + b + 3)(b + 3) + 9A 2 )t + C 2 ) 2 ' 


(3.11) 


where we neglected Ci and used equations (2.5), (3.9) and (3.10). So, using the equation 
(2.7) we can investigate evolution of the deceleration parameter. At the very early universe 
one can obtain the following expression, 


q = -l+A-Bt 2 + 0(t 3 ), 


(3.12) 
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where constants A and B defined as 


A = 
B = 


3A + b + 3 
2 ’ 

3b 2 (3 A + 6 + 3) 3 


(3.13) 


It is clear that — 1 + A is positive for positive values of b and A, hence the value of q is 
positive at the early universe which is expected. 


3.2.2 The case n = 2 

Again, we consider the early universe cosmology where energy density is very high and 
assume that a = 0.5 and n = 2. Under these assumptions the ECG behaves as barotropic 
fluid with quadratic EoS. In that case we can obtain analytical solutions of the equations 
(2.8) and (2.9). First, we can obtain the following expression from the equation (2.8), 


3\/2Atan 


-i 


t — 


3Ap eC g \ 
3A+6+3 I 


Ve 


(3 A + 6 + 3)2 


(3A + b + 3 )yj Pecg 


+ C3 — 0, 


(3.14) 


where C 3 is an integration constant. At the initial time with the high density we can 
obtain, 

3A + b + 3 , 2 


Pecg 


3A 


■ tan t, 


(3.15) 


where t = v ^( 3y4 +^+ 3 ) 2 ( t+Ca I i s use d. In that case one can obtain, 


P<f> 


\/6b(3A + b + 3) (t — tan t) + C 4 
324A 2 ’ 


(3.16) 


where C 4 is another integration constant. 

Easily we can check that the Hubble expansion parameter is decreasing function at the 
initial times. Also, analysis of the deceleration parameter tells that the universe may begin 
with the negative q, if value of b be smaller than one, then transfer to the universe with 
positive q and again return to the deceleration phase with the negative q. On the other 
hand, with b > 1 we have initial universe in acceleration phase which transferred to the 
deceleration phase at the late time. 

Behavior of the interaction term Q is interesting, it has periodic like feature during time. 
It begin with a high value at initial time and decreased to a small value and again grows 
up to a high value, and repeat this procedure. However, it is corresponding to the early 
universe. 

It is useful to investigate behavior of the scalar field. We can use some assumption about 
the scalar potential V(4>) from the previous works such as the PNGB model [93] with a 
phantom potential given by, 


V(0)=pw-M\l 



(3.17) 
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where M and / are constants and o is present value of density which can considered as a 
constant. There is also other model with a local phantom potential including the Gaussian 
potential [37], 

V{(j)) = pw-M\l-eb), (3.18) 

where a is a constant. Moreover the quadratic phantom potential may be written as [37], 

V(4>) = PP o - V 2 4> 2 , (3.19) 

where V 2 is a constant. All of the above models for the infinitesimal values of the scalar 
fields behave as the following, 

G((/> <C 1) = ai + a 2 cj) 2 + 0((j) 3 ), (3.20) 

where a\ and a 2 are arbitrary constants. In that case we can obtain typical behavior of 
the scalar field in terms of time corresponding to our model as illustrated in Fig. 2. We 
can see that the scalar field is increasing function of time, and its evolution in the early 
universe is faster than that in the late time. 



Figure 2. Scalar field <f> in terms of t, for A = 1/3, C 3 = 1, C 4 = 1, a = 1/2, b = 1, ujq = — 1 and 
n = 2. 

3.2.3 The case n > 2 

For the larger n, at high density regions, we can consider only the last term of the expansion 
in EoS given by the Eq. (2.4). In that case the solution of the equations (2.8) and (2.9) 
obtained easily as follow, 

2 

2n — 1 

, (3.21) 


Pecg — 


y/6A(2n - 1 )t + 2 C 5 
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and 


P<t> 


_V6A(2n - l)t + 2C 5 


x 


A 2 b 2 (2n — l ) 2 2 

(A(2n - 3) + C 6 ) 2t ’ 


(3.22) 


where C 5 and Cq are integration constants and should be considered as C 5 <C 1 and C& « 1 
to satisfy p ecg S> p$ condition, so we can neglect C5. 

We find that the Hubble expansion parameter is decreasing function of time which decreased 
also by increasing n. Also, we can investigate evolution of the deceleration parameter and 
see that it begin with positive value and yields to -1 at the late time (see Fig. 3), which 
shows deceleration to acceleration phase transition. 



t 


Figure 3. Deceleration parameter q in terms of t for the case of larger n. 


3.3 Special solution 

In the previous subsections we restricted ourself on the case of higher ECG energy density 
than scalar field energy density. This assumption simplified our equations and help us to 
obtain some analytical expressions. Now, we try to study more general case without any 
condition on the energy densities. In that case, from the conservation energy equations, 
one can obtain, 

Pecg + 3Ff(l + Uj ecg )pecg + 3 bH S = 0, (3.23) 

and, 

P</) + 3-ff(l + uj<f>)p<p — 3bH 3 = 0. (3.24) 

We can consider above relations as cubic equation for H and solve them to obtain, 


H = 


X 3 2(1+ UJecg ) Pecg 

6 b x 5 

Y 5 2(1 +u < j > )p ( f > 

/"* 7 I 1 5 


(3.25) 
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where we defined 


X = 

Y = 


-36 p ecg + 12y/3\ 


3bp 2 ecg + 4/?| C3 ( 1 + 3u} ecg (l + UJecg) + UJ^ cg ) . „ 


^ + 12 +»<,)+*$> 


(3.26) 


Therefore, having /3 ecg , can gives us Hubble parameter as well as p eC g- Then, from the 
equation (2.5) one can obtain p^ and use it in the equation (2.6) to find p^. Moreover, we 
can use the equations (2.3) to obtain scalar field and scalar potential. So, we can use some 
ansatz for p ecg to have all cosmological parameter and compare them with observational 
data. 

Here, we assume that the extended Chaplygin gas has the following form, 


Pecg 



(3.27) 


where a 2 and m are free parameters which will be fixed using recent observational data. 
The ansatz (3.27) tells that p ecg —>• 00 at the early universe and p ecg —>• y/a\ at the late 
time. These are agree with the fact that the energy density should be decreasing function 
of time in an expanded universe. We expected that the parameter a± be small to have 
small energy density at the late time. This yields to a Hubble parameter (as well as (f>) 
as decreasing function of time which is expected. Also, we can find that the deceleration 
parameter yields to -1 at the late time, with possibility to have acceleration to deceleration 
phase transition. All of them depend on choosing appropriate parameters. We will fix 
them using recent observational data in section 4. Interestingly, we can analyze, 

n ecg = a 2 , (3.28) 

and, 

= ^ 2 fl2 ’ ( 3 ‘ 29 ) 
to see the exact equality H ecg = It is illustrated in the Fig. 4. For various n the curves 
of H ecg and have overlap and recovers each other exactly. It will be happen even for 
n = 1, therefore interaction term may be cause of this event. 


3.4 Inflationary solutions 

Finding H approximately as a constant, we can investigate inflationary solutions in the 
model of interacting extended Chaplygin gas with phantom field. Recently, an inflationary 
model of modified Chaplygin gas with quintessence field has been studied [94] . The present 
work may be extension of that paper to the case of ECG in the context of phantom 
cosmology. Therefore, there are three main differences between this work and Ref. [94]: 
The first: extension of MCG to ECG. 

The second: exchange of quintessence field to phantom field. 
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Figure 4. Typical behavior of Q ecg and in terms of cosmic time, n = 1 blue, n = 2 green, 
n = 3 yellow, n = 4 red. 


The third: consideration of interaction term. 

In order to obtain such solutions we combine equations (2.8), (2.9) and (2.5) to find the 
following differential equation, 


H -ih + 2 H ^ +b) - P i 


+ 




LH 2 - **]” * 


B 


HH 2 -ptY 


= 0 . 


(3.30) 


The simplest inflationary solutions obtained under assumption of constant p^. Fig. 5 
shows that Hubble parameter obtained as a constant if we assume p^ = C , where C is an 
arbitrary constant. We exam our model for n = 1, 2,3 and found that the Hubble expansion 
parameter is a constant during cosmic time. Therefore, We can find ECG energy density 
as the following constant, 


Pecg — C. 


(3.31) 


3.5 Bouncing solutions 

An interesting solution of the singularity problem of the standard cosmology known as 
bouncing solution [89, 95-100]. The scale factor is decreasing in the contracting phase, 
and increasing in the expanding phase. It means that H < 0 transits to H > 0. We can 
reconstruct scalar potential as logarithmic function with behavior represented in the left 
panel of the Fig. 6. Then, we can obtain behavior of the scalar field with time (see right 
panel of the Fig. 6). 

Therefore, we can obtain behavior of the scalar potential in terms of the scalar field in 
the Fig. 7. We find that the scalar potential is combination of some logarithmic functions 
and totally it is decreasing function of scalar field. 

The scalar potential has complicated and long expression with logarithmic nature. 
Motivated by this behavior we can reconstruct scalar potential with the simple logarithmic 
form, 

V{4>) = d\ In [d 2 (0 — 0o)L (3-32) 
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Figure 5. Hubble parameter versus time for a = 1, B = 3.4, A n = ^ and C = 85. n = 1 (solid 
line), n = 2 (dashed line), n = 3 (dotted line). 




Figure 6. Typical behavior of scalar potential and scalar field in terms of cosmic time for a = 0 . 5 . 

where d\ and ^2 are arbitrary constants. This is very simplified potential but gives the 
similar results. We can use it to obtain p^ and therefore, <fi. Then, we can find Hubble 
parameter via the relation, 

H = (3.33) 

where V' = It yields to the Hubble parameter represented by the Fig. 8. We can see 
that the Hubble parameter is negative for t < 1 with u)<f, > — 1 and it is positive for t > 1 
with < — 1 (see Fig. 8). 

It yields us to the following general solution for the ECG energy density, 

Pecg = h \nb 2 t, (3.34) 

where b\ and 62 are constant parameters related with d\ and c? 2 - In that case, using the 
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Figure 7. Typical behavior of scalar potential in terms of scalar field for a = 0.5. 



Figure 8. Hubble parameter versus time for d± = 5 and d ,2 = 0.965. 

equations (2.8), (2.9) we can obtain behavior of the Hubble expansion parameter (see Fig. 
9). It means that by using logarithmic behavior of the scalar potential and ECG energy 
density the scale factor decreased first, corresponding to negative value of the Hubble ex¬ 
pansion parameter, and increased after some time corresponding to positive value of the 
Hubble expansion parameter. 

We can see that in the cases of n = 1,2,3, the Hubble parameter is completely in¬ 
creasing function of time which yield to a constant at the late time. The transition form 
H < 0 to H > 0 happen at t = 1 for selected values of parameters. In the case of n = 4 the 
Hubble parameter decreased first to a minimum negative value, then increases to positive 
value. The phase transition time is similar to the previous cases. Cases with higher n 
yield to initially positive H which suddenly take negative value after the early universe 
and behave as previous cases after the initial time. It is obvious that, at the late time there 
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Figure 9. Hubble parameter versus time for a = 1, B = 3.4, A n = 4, b\ = 1 and 62 = 1. n = 1 
(solid blue line), n = 2 (dashed green line), n = 3 (dotted yellow line), n = 4 (dash dotted red line), 
n = 5 (solid cyan line). 


is no differences between various n. The fact is that, importance of the first term in the 
right hand side of the equation (2.4) is in the early universe. After some time, the second 
term is dominant. 

In that case the scalar field energy density is given by, 


Ptf, = 3H 2 - 61 In b 2 t, 


(3.35) 


where, 


H = -b\ \nb 2 t 
3 


1 +nA{b 1 In b 2 t) n ~ l - 


B 


( b\ \nb2t) l+a 


-1 + 


1 - 


\ 61 In (M ) 2 ((1 + nAb\ In b 2 t) n 1 - (6i ln g t) i+ a ) 

and we assumed A„ = A. 


(3.36) 


4 Perturbations 

4.1 Scalar perturbations 

In this section we give density perturbation analysis of our model. Already, density per¬ 
turbation of a universe dominated by Chaplygin gas was studied by the Ref. [101]. Now, 
we use their results to write the following perturbation equation, 

3 C 2 

5 + H[2- 3(2w - C 2 s )]6 - -H 2 (l - 6 C 2 - 3w 2 + 8 c j)5 = (4.1) 
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where 5 is a density fluctuation, k is the wave number of the Fourier mode of the pertur¬ 
bation, and 

C s 2 = (4.2) 

P 

is squared sound speed. 

Analytical solutions help us to investigate density perturbation and stability analysis of 
the model. We can investigate stability of the model by using squared sound speed. If 
Cg > 0 then, the model is stable. We sould note that it is different with the gravitational 
instabilities allow for structure formation as observed in the sky. 

In the simplest case, n = 1 we can find that our model is stable at the initial time (see Fig. 
10). We can see from Fig. 10 that the case of non-interacting (b = 0) yields to constant 
sound speed. By increasing strength of interaction, the sound speed decreased and yields 
to a constant at the late time. So, the model is completely stable for n = 1. 




t t 


Figure 10. Squared sound speed in terms of t for A = 1 / 3 , C\ -C 1, C 2 = 1, a = 1/2 and n = 1 . 
Left: 5 = 1 , = 0.5 (dash dotted line), = 0 (dotted line), = — 0.5 (dashed line), = —1 

(solid line). Right: w = — 1 , b = 0.5 (dash dotted line), 5 = 0 (dotted line), b = 0.2 (dashed line), 
5=1 (solid line). 

In the case of n = 2 of the early universe we can obtain stability condition. Fig. 11 
shows that there are some instabilities depend on value of the interaction strength. How¬ 
ever, instabilities vanished soon and our model will be stable. Opposite of the previous 
model with n = 1, the case of b = 0 in present model with n = 2 don’t lead to the constant 
sound speed. 

For the case of n > 2 we can see behavior of squared sound speed in Fig. 12. We can 
see that presence of interaction makes some instabilities at the initial times. So, in the 
case of b = 0 the model is completely stable at the early universe. However, for the small 
values of b (b -C 1) we have stable model at the initial stage. 

Numerically, we can obtain behavior of 8 for the cases which studied in this paper. In 
the simplest case of n = 1 we see behavior of 8 in the Fig. 13. It shows that perturbations 
grow initially and vanished at the early universe. There is a maximum value at the early 
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Figure 11. Squared sound speed in terms of t for A = 1/3, C 3 = 1, C 4 = 1, a = 1/2, = — 1 

and n = 2 . 6 = 2 (dotted line), b = 1.2 (dashed line), 6=1 (solid line), 6 = 0 (dash dotted line). 




Figure 12. Squared sound speed in terms of t for A = 1/3, C 5 -C 1, C 6 = 1, a = 1/2 and w = —1 
in the case of n > 2. Left: n = 10, 6 = 0 (dash dotted line), 6 = 0.4 (dotted line), 6 = 0.6 (dashed 
line), 6 = 0.8 (long dashed line), 6=1 (solid line). Right: 6 = 0, n = 3 (dash dotted line), n = 6 
(dotted line), n = 8 (dashed line), n = 10 (solid line). 


universe. 

Fig. 14 shows behavior of 5 corresponding to n = 2 at the early universe under assumption 
of Pecg p<t>- It is clear that the perturbations vanished suddenly at the initial times. 
The maximum value is at the initial time exactly and perturbation is totally decreasing 
function of cosmic time. 

For the case of n > 2 we can see from the Fig. 15 that value of interaction coefficient is 
important in evolution of the perturbations. In the case of b = 1, value of 5 reaches to the 
large number, however vanished at the early universe. In the case of without interaction 
(6 = 0 ), perturbation 5 is decreasing function of time and vanished rapidly. 

Therefore, we show that, density perturbations vanished at the early universe and we can 
have stable model at the late time. 
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Our numerical results are base on the simplest cases with k = 0, however we can obtain 
similar results for the cases of k 7 ^ 0 . 



Figure 13. Time evolution of 6 for the simplest case: p ecg pp and n = 1 at the early universe. 



Figure 14. Time evolution of <5 for the first exited case: p ecg pp and n = 2 at the early universe. 


4.2 Tesor-to-scalar ratio 


We can use our solutions to investigate the tensor to scalar ratio and compare it with recent 
observational data [102], We follow the Ref. [103] to write, 


r « 32 



(4.3) 


- 17 - 






(a) (b) 




Figure 15. Time evolution of 5 for the case of large n and p ecg p$ at the early universe, (a) 
6 = 0 and 6 = 2 represented by solid line and dashed line respectively, (b) the case of 6 = 1. 

We use energy density given by the equation (3.27) to find behavior of the tensor to scalar 
ratio. Results of our numerical analysis represented in the plots of the Fig. 16. For 
simplicity we take B = a\ = 02 = A n = 1 and a = 0.5, hence the only free parameter is m 
(see Eq. (3.27)). We can see that higher n are closest to the observational data. 

We can find approximately similar behavior for bouncing solution given by ansatz (3.34). 
In this case also having agreement with KBP observation [102] and obtaining r < 0.09 
needs to include higher n. 



t 


Figure 16. The tensor to scalar ratio in terms of cosmic time with B = A n = 1, ai = a 2 = 10 , 
6 = m = 0 . 1 , and a = 0 . 5 . n = 1 blue, n = 2 green, n = 3 yellow, n = 4 black, n = 5 red. 

We can find approximately similar behavior for bouncing solution given by ansatz 
(3.34). In this case also having better agreement with r < 0.09 condition for higher n. 

5 Conclusion 

In this paper, we proposed a cosmological model based on phantom dark energy which in¬ 
teracts with the extended Chaplygin gas. Therefore, we have two-component fluid universe 
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with possibility of interaction between components. First of all we introduced our model 
and wrote basic equations. Then we discussed about solutions for the several cases and 
found inflationary and bouncing solutions. Numerically, we found behavior of the Hubble 
parameter at the late time for the general case and discussed briefly about the special case 
of n = 1 which recovered modified Chaplygin gas. At the early universe we discussed cases 
of n = 2, n = 2 and n > 2 separately and found agreement with the expected behavior 
of cosmological parameters. In the case of n = 2 we reconstructed phantom potential and 
obtained behavior of phantom field as a function of cosmic time. 

We also analyzed density perturbation of the model. By using squared sound speed we 
investigated stability of the model. In the case of n = 1 we found that our model is stable 
at the initial time. We have shown that the case of non-interacting yields to constant sound 
speed. By increasing strength of interaction, the sound speed decreased and yields to a 
constant at the late time. We found some instabilities at the early universe of the case 
n = 2 depend on value of the interaction strength. However, instabilities vanished soon 
and our model will be stable. For the general case of n > 2, we have seen that presence of 
interaction makes some instabilities at the initial times. However, for the small values of 
interaction coefficient we have stable model at the initial stage. Finally, comparing with 
recent observations of KBP [102] suggests that presence of higher order terms of the ex¬ 
tended Chaplygin gas are necessary to obtain agreement with observations. 
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